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We study the spin- 1/2 trellis lattice Heisenberg model, a coupled spin ladder system, both 
by perturbation around the dimer limit and by quantum Monte Carlo simulations. We discuss 
the influence ol the inter-ladder coupling on the spin gap and the dispersion, and present re- 
sults for the temperature dependence of the uniform susceptibility. The latter was found to be 
parameterized well by a mean-held type scaling ansatz. Finally we discuss fits of experimental 
measurements on SrCu2C>3 and CaX^Os to our results. 
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§1. Introduction 

Spin- 1/2 ladder models 1 ) can describe the magnetic 
behavior of a variety of quasi-one dimensional materi- 
als. Examples include the cuprate materials SrCu203 2 ) 
and LaCu02.5 3 ' ) and the vanadate CaV205. 4 ) Spin ex- 
citations in isolated ladders have a finite energy gap, 
which makes them prototypical spin liquids. 1 -' This is of 
interest in relation with high temperature superconduc- 
tors, since upon doping they become doped resonating- 
valence-bond liquids, with a spin excitation gap and dom- 
inating quasi-long range pairing correlations. 

While isolated ladders are relatively easy to study and 
offer a variety of interesting phenomena, it is still nec- 
essary to study the inter-ladder coupling. Although the 
inter-ladder coupling is weak it is still necessary to ob- 
tain superconducting long range order, as observed in the 
ladder compound (Sr,Ca)i4Cu2404i. 5 ) Even in undoped 
ladder systems the weak inter-ladder coupling can play 
an important role. In LaCu02.5 & n inter-ladder coupling 
of only one tenth of the intra-ladder coupling is sufficient 
to destroy the spin gap of the isolated ladder and leads 
to antiferromagnetic long range order. 6 ) 

One of the authors fitted the uniform magnetic suscep- 
tibility of SrCu203 to that of an isolated ladder model 7 ) 
and found that the best fit is achieved with a ratio J±/J\\ R 
0.5, where Jj_ is the coupling along a rung of the ladder 
and Jy that along the legs of the ladder. This result 
is surprising since the Cu-Cu distance across a rung is 
shorter than that along a leg and the rung coupling J± 
thus expected to be larger than the leg coupling. It was 
suggested that this apparent ratio of J±/J\\ ~ 0.5 may be 
due to the neglect of inter-ladder interactions, although a 
mean-field analysis of the influence of these interactions 
suggested a similar ratio. In this paper we report on 
an investigation of the effect of the inter-ladder coupling 
on the magnetic susceptibility and the magnon disper- 



sion of coupled ladders and dimers. Fits of the results to 
experimental measurements reveal that the inter-ladder 
coupling does not change the original estimates signifi- 
cantly and will be discussed in detail in a forthcoming 
publication. 8 ) 

The susceptibility of another ladder-like compound CaV205 
was fitted by Onoda and Nishiguchi to an isolated dimer 
model with Jy = 0. 9 ) Again this result is surprising as 
the bond lengths are similar. We perform simulations 
also for this compound to estimate the strength of inter- 
dimer couplings. 
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Fig. 1. The lattice structure of the trellis lattice Heisenberg 
model. J± is the exchange constant across the rung, Ja is that 
along the legs and J' is that between ladders. 



Figure 1 shows the structure of the trellis lattice. It 
consists of spin ladders (with coupling along the legs Jn 
and across the rungs Jj_, coupled by frustrating inter- 
ladder "zig-zag" couplings J'. The size of the lattice 
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used in the calculations is denoted by N x x N y , where 
N x is the number of the spin along the x-axis and N y 
along the y-axis. 

In SrCu2C>3 and CaV20s the interactions in the ladder 
( Jl and Ju ) are antifcrromagnetic. The inter-ladder cou- 
pling (J') in SrCu2C>3 is expected to be ferromagnetic, 
since it occurs via 90° Cu-O-Cu bonds. In CaV20s the 
local chemistry is more complex, and the sign of the ex- 
change coupling J' has not been definitively established. 

§2. Spin Gap and Magnon Dispersion 

In this section, we discuss the influence of the inter- 
ladder coupling J' on the spin gap and the magnon dis- 
persion. Following refs., 10, but extending their ideas 
to include the inter-ladder couplings, we start from the 
coupled dimer system, where J_l » Ju, J'. 

The ground state in the dimer limit (J|| = J' = 0) 
consists of spin singlets on the rungs, with an energy per 
rung of E S i n gi e t = — | Jj_. Including the inter-dimcr cou- 
plings perturbatively, the ground state energy per rung 
to second order is 

f 3 r 3 J| i 3J ' 2 ran 

The lowest lying excitations arc spin-1 magnons, where 
one of the rung singlets is turned into a triplet. These 
magnons disperse due to inter-dimer couplings. The low- 
est branch is 
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uj(k±, fen 



i + 4Ji 8 Jl 



1 J 

+J||Cos(fc||) - - JLcos(2fc||) 



+ HJ'\- 



J' 



\JV h\J' 



2J_l 



2J± )cos( T )cos( y ) 



— — (cos(fcj_) cos(fc||) + cos(fc^) + COs(fc||)) 



+■ 



8J 
Jul J' 



2J ± C0S ( T )C ° S( — } ' 



(2.2) 



where we have chosen units such that the lattice param- 
eters (see Fig. 1) are a = b = 1. 

Minimizing uj(k±, k\\), we obtain for the the spin gap 
A s in next to leading order 




J||>|0.25J'| 
IJ'I Ju < |0.25J'| 



The minimum of the dispersion is at momenta 

fc_L = 
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We can extend above results to arbitrary ratios of 
J\\/J± by replacing the second order expression of the 
dispersion due to Ju 
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by the exact dispersion e{k\\). 

Substituting the J|| terms in eq. (2.2) by e(fc||) we ob- 
tain 
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+ cos(fc||)]. (2.7) 
Expanding e(fc||) around the minimum at fc|| = ir as 

(2.8) 



e(7r + fc||) = ^A 2 +v 2 k ll 2 , 



where Ao is the spin gap for the spin ladder system and 
v is the magnon velocity we get a minimum in the dis- 
persion at 

A Ap 2 J' 



(2- 



Jj 
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and for the influence of the inter-ladder coupling J' on 
the spin gap 



where 



C 



(2.10) 



(2.11) 



The change in the spin gap is small, second order in J' 
and with a typically small prefactor |(— ) ~ 0.1. 10 ' 12 ) 
We have calculated the spin gap on (0, tt) on 4 x 4 and 
4x6 lattices by exact diagonalization using the Lanczos 
algorithm 13 ) and observed only minimal changes in the 
spin gap, in agreement with above arguments. These 
system sizes are however too small to allow quantitative 
comparisons. 

§3. The Uniform Susceptibility 

3.1 Weakly Coupled Ladders 

While the inter-ladder coupling J' (J' < 0.2Jm) has 
negligible influence on the spin gap it still modifies the 
uniform magnetic susceptibility \ at intermediate and 
high temperatures. Using the quantum Monte Carlo 
loop algorithm 14, 15 ) with improved estimators 16 ) we cal- 
culated the temperature dependence of \. 

We considered ladder systems with both J±/J\\ = 1 
and the experimentally relevant coupling Jj_/J\\ = O.5. 7 ) 
The inter ladder couplings were chosen to be J' / J\\ = 
±0.1, ±0.2, ±0.5 for the isotropic ladders and J' / J\\ = 
±0.1 and ±0.2 for Jj_/J|| = 0.5. 

The quantum Monte Carlo simulations suffer from a 
negative sign problem due to the frustrated inter-ladder 
couplings J'. This sign problem was alleviated a bit by 



Quantum Monte Carlo Simulation of the Trellis Lattice. 



3 



using improved estimators both for \ an d f° r the av- a fourth-order high temperature expansion 



erage sign. 16 ) Still the sign problem restricts the QMC 
simulations to rather small lattices and not too low tem- 
peratures. In the temperature regime where QMC simu- 
lations were feasible a lattice size of 4 x 16 was sufficiently 
large and our results are not biased by finite size effects. 
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Fig. 2. Temperature dependence of the uniform susceptibility \ 
of the trellis lattice Heisenberg model with J±/J\\ = 1 and 
J' /J\\ = 0, ±0.2 and ±0.5. The dotted lines are fourth-order 
high temperature expansions. 
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Fig. 3. Temperature dependence of the uniform susceptibility \ 
of the trellis lattice Heisenberg model with Jx/J|| = 0.5 and 
J 1 /J\\ = 0, ±0.2. The dotted lines are fourth-order high temper- 
ature expansions. 
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Fitting magnetic susceptibility measurements on spin 
ladder materials is much simplified if approximate ana- 
lytic expressions for the magnetic susceptibility are avail- 
able. The single ladder susceptibilities have been param- 
eterized before. 7, 8 ' 12, 

To extend these parameterizations to include the inter- 
ladder coupling we follow a mean field-type scaling ansatz 
(MFTS) of ref. 7 ) and scale the trellis lattice susceptibility 
to the susceptibility xi(J±/J\\) of a single ladder: 



with 



J \\ J \\ 



j' 

!( T ) 
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J\\ 



(3.2) 



(3.3) 



and A and B as temperature independent constants. At 
low temperatures this ansatz gives the same gapped be- 
havior as a single ladder, which is reasonable since the 
gap is not changed much. Fixing A = 2 also recovers the 
correct high temperature Curie- Weiss law. By fitting our 
QMC results in the temperature range 1.0 <T/Jj_ < 1.5 
to eq. (3.2) we obtain 



B = 0.3436(3). 



(3.4) 



As a check of the quality of the MFTS ansatz we com- 
pare, in Fig. 4, the single ladder susceptibilities obtained 
from the QMC data of coupled ladders by inverting eq. (3.2) 
. We find that the ansatz is useful for the whole tem- 
perature range as long as the inter-ladder coupling J' is 
small. 

3.2 Coupled Dimers 

We can repeat a similar analysis for weakly coupled 
dimers instead of weakly coupled ladders. We considered 
leg couplings of J\\/J± = 0, 0.1 and 0.2 and inter-ladder 
couplings J'/J± = 0, ±0.1 and ±0.2 on systems of size 
32 x 16, where finite size effects are again negligible. In 
Fig. 5 we show the uniform susceptibility for J\\/J± = 
0.1. 

We again make use of a MFTS ansatz to obtain an 
approximate analytic expression for the susceptibility 



(3.5) 



Our results are shown in Figs. 2 and 3 together with 
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Fig. 4. (a) Scaling plot of weakly coupled ladders, (b) differences 
between the scaled and actual susceptibilities. 



with 



~ V ~ V ~ V 



(3.6) 



and A and B temperature independent constants. The 
high temperature expansion again fixes A = 2, and fits 
to the QMC susceptibility data in the temperature range 
1.0 < T I Jj_ < 1.5 give 



B = 0.6940(3). 



(3.7) 



Scaling plots (Figure 6) show that the agreement is 
not as good as in the coupled ladder case. This is not 
very surprising, since, as can be seen in eq.( 2.3) and 
cq.( 2.10), the spin gap of weakly coupled dimers de- 
pends more on the inter-dimer couplings than that of 



Jh/Ji =0.1 




Fig. 5. Temperature dependence of the magnetic susceptibility 
of the trellis lattice Heisenberg model in the regime of weakly 
coupled dimer. The coupling ratios are J\\/J±_ = 0.1 and 
J'/J± = ±0.1 and ±0.2. The dotted lines arc again the fourth 
order high temperature expansions. 



weakly coupled ladders. Still the MFTS ansatz gives a 
reasonable analytic parameterization of the susceptibili- 
ties of weakly coupled dimers. 

§4. Comparison with Experiments 

In this section we wish to briefly compare our QMC 
results with experimental measurements of the suscepti- 
bility of SrCu203 and CaV205. A forthcoming publica- 
tion 8 ) will present detailed fits and comparisons. 

4.1 SrCu 2 3 

With regard to fits of the experimental susceptibility 
measurements on SrCu203 our main result is that the 
inclusion of inter-ladder couplings does not modify the 
previous estimate J±/J\\ ~ 0.5 7 ) considerably. We can- 
not determine the value of J' from these fits since in the 
experimentally accessible temperature range T < 650K 
the dependence of the susceptibility on J' is weak. 

4.2 CaV 2 5 

It was proposed previously that the magnetic suscepti- 
bility of CaV2C>5 can be fit by that of dimers. 9 ) We have 
performed a series of fits on new experimental data of 
Isobe and Ueda. 18 ) 

One problem is that the samples are not pure CaV205 
but contain a few percent CaVaOy. By X-ray struc- 
tural analysis Isobe and Ueda determined that the cur- 
rent sample contains 4.1% CdN^O^. By subtracting the 
separately measured susceptibility of CaV307 18 ) we ob- 
tained experimental data for pure CaV 2 05. 

The main result of our fits is that the inter-dimer cou- 
plings J|| and J' are both much smaller than Jj_, and we 
can give some constraints. 
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To determine the spin gap A we fitted the low tem- 
perature experimental data to the expression for a spin 
ladder at T < A 12 ) 

x( T ) = ^Q+Xo + ^cxp(-A/T), (4.1) 

where the g-factor of Vanadium was determined to be 
g = 1.96 by ESR measurements. 9,18 ^ The fit parame- 
ters were C = 1.6 x 1CT 3 cm 3 K/mol V, 6 = -8K and 
Xo = 1-23 x 1CT 5 cm 3 /mol V. We use this fit to sub- 
tract the first term, which is a Curie- Weiss term due to 
paramagnetic impurities, and the second, temperature 
independent, constant contribution. 

Next we fit the low temperature data (T < 200K) to 
QMC results for an isolated ladder for J\\/J± — ±0.1, 
±0.2 and ±0.3. This fit is reasonable because we will see 
later, that we are in the ladder case \J'\ <C J\\, where the 
influence of J' on the low temperature susceptibility is 
negligible. In all cases we found Jj_ w 670K. 




100 200 300 400 500 600 700 
T(K) 



Fig. 7. Fit of the uniform magnetic susceptibility for CaV^Os, 
assuming J\\/J± = 0.1. The fit parameters are J±_ = 672K, 
J|| = 67K and J' = 45K 

Next, for each of the above ratios of J|| / J± with J'/ J± = 
we use the MFTS eq. (3.5) to fit the susceptibility in 
the range 200K < T < 700K and determine J'. Good 
fits are obtained only for Ju/J± =0.1 and 0.2, but not 
for J\\/J± < or J\\/J± = 0.3. Assuming J\\/J± = 0.1 
we obtained J± w 670K, J {1 w 67K and J' w 45K. This 
fit is shown in Fig. 7. For the other ratio J\\/J± = 0.2 
we got an equally good fit with J± w 665K, J|| w 135K 
and a ferromagnetic J' w — 25K. 

Given the good quality of both fits it is hard to deter- 
mine the exact values of the couplings Jn and J' from fits 
to the uniform susceptibility alone. But we can give some 
estimates and constraints: J± ~ 670K, OK < Jm < 200K 
and J' + J|j w 110K. The dispersion relation of the 
magnons however depends sensitively on the ratio J|| / J' 
and could offer a way to determine these couplings more 



precisely. 

§5. Conclusions 

We have studied the trellis lattice Heisenberg model 
by quantum Monte Carlo and a perturbation expan- 
sion around the dimcr limit. We confirmed that for 
weakly coupled ladders the influence of the frustrated 
inter-ladder coupling on the spin gap is small. 

We calculated the uniform susceptibility by QMC sim- 
ulations and found that a mean field-type scaling ansatz 
gives a reasonable analytic parameterization of its tem- 
perature dependence in the whole temperature range. 

This parameterization was in turn be used to fit exper- 
imental measurements on the compounds SrCu203 and 
CaV205. For SrCu203 it was found that the inter-ladder 
coupling does not significantly modify the previous es- 
timate of the intra-ladder coupling ratio J±/Ju ~ O.5. 7 -* 
CaV205 was confirmed to be a weakly coupled dimer sys- 
tem with inter-dimer couplings about an order of mag- 
nitude smaller than the intra-dimer coupling. 
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Fig. 6. Scaling plots for the weakly coupled dimer regime for couplings (a) J\\/J±_ = 0, (b) Jy/Jx =0.1 and (c) J\\/J±_ = 0.2. (d) 
shows the differences between the scaling curves 



